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Abstract 

Several simple approximations for the evolution during the nucle- 
ation period have been presented. All of them have been compared 
with precise numerical solution and the errors have been estimated. 
All relative errors are enough small. 

In the process of the first order phase transition which is ordinary consid- 
ered on the example of condensation, i.e. the formation of a liquid droplets 
phase in a mother metastable vapor phase, a period of embryos formation 
plays the main role. This period is called a nucleation period. Despite the 
variety of substances one can extract some common features in evolution 
during this period. These features form the base of the theoretical descrip- 
tion. The most precise solution of this problem was given in [T| and is based 
on the universality of the nucleation period. But this universality disappears 
when we consider more specific phenomena and more complicated situations. 
That's why we need some simple approximations which will allow us to in- 
vestigate more complicate situations. The goal of this publication is to give 
these approximations. Beside the simple functional dependencies we shall see 
the physical reasons which lie in the base of these approximations. This gives 
us an opportunity to give more simple model for the evolution during the 
nucleation period which can be used in consideration of more specific char- 
acteristic of the phase transition, for example, for manifestation of stochastic 
effects of embryos appearance. 
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1 Main equations 



As an example we consider the case of free molecular regime of growth of 
embryos in a three dimensional space. Then the rate of growth of the embryo 
with v molecules inside in time t is given by a standard relation 



Here r is some characteristic time and we extract the dependence on the 
density of the vapor phase via supersaturation ( which is 



where n is the vapor molecules density, is the same characteristic for the 
saturated vapor. 

The value of supersaturation ( characterizes the power of metastability 
in the system. The analogous value which would be attained in the system 
without embryos formation will be denoted by This value is completely 
determined by external conditions and, thus, is supposed to be known. 

The distribution of the embryos / over their sizes can be taken as the 
stationary one and the following approximation for / dependence on ( can 
be written [I] 



Here and later index * marks values at some characteristic moment £* which 
can be chosen as the moment when the half of the total number of droplets 
is formed. The value T is the derivative of the critical embryo free energy 
over the supersaturation multiplied by $*. 

During the nucleation period one can use linearization of the ideal super- 
saturation and after the appropriate choice of scale 



Here z is rescaled time (see pQ). One can introduce x ~ z — ap dr (a 
is a constant 0) and establish correspondence between z and x. Under 
the collective regime of vapor consumption which takes place under the free 
molecular regime the conservation law for a substance can be written as |T] 



dt 



C/r, p dr = v 



.1/3 



/(C) = /(*.) ex P (r 



$ = + z 
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with the appropriate value of the constant c = 0.189. Then the curve exp(x — 
g(x)) has maximum at x = (this can be regarded as a choice of c). 

Now we shall construct approximations for solution / ~ exp(x — g(x)) of 
(P). It is shown in Figure 1. 
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Figure 1: Precise solution for the size spectrum 



This curve has a universal form which has been noticed in P . 

2 Iteration approximations 

The iteration solution of eq. was given in 1 j^j. The recurrent procedure 
is defined as 

gi + i — c (z — x) 3 exp(x — gi{x))dx 

J — oo 

with zero approximation go = 0. The first iteration for f(x) ~ exp(x — g(x)) 
has the following form 

fx = exp(x — pcexp(x)) 

1 Solution in [2] requires the new set of characteristic parameters at every step of iter- 
ation procedure. 
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p — x 3 exp(— x)dx = 6 
Jo 

pc = 1.134 

One can see that the first iteration is very close to the real precise solution 
which can be seen from Figure 2. The upper curve is the precise solution, 
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Figure 2: Precise solution and the first iteration for the size spectrum 

the lower one is the first iteration. 

One can see that the first iteration is already a good approximation. Now 
we shall formulate some new approximations with more simple physical sense. 

3 Similarity of the embryos formation condi- 
tions 

Under the conditions of the metastable phase decay one can see in jl] the 
similarity of conditions under which the embryos are appearing in the system. 
It was the consequence of the scale transformation of time. Here the situation 
seems to be another. The time explicitly stands in the term exp(x) for 
/ ~ exp(x — g(x)) which characterizes the action of external forces. Let 
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us consider the subintegral expression in the integral calculated for the first 
iteration. The size pd r of the embryo can be expressed as 

Pdr = T^Idt {z ~ X) 

To observe the similarity we need to see the subintegral equation in the 
variable p = z — x. The subintegral function g su b i in the expression for g\ is 
the following 

9sub 1 ~ exp(2)p 3 exp(-p) 

Factor exp(z) is not more than a scaling factor and we have the universal 
expression 

g S ub 1 ~ p 3 exp(-p) 

for subintegral function in the first iteration. Because the first iteration 
resembles the precise solution we see that the last relation means that condi- 
tions for appearance of the droplets in different moments of time are similar, 
only the scaling factors change. Here it isn't clear how to use this property, 
but in the direct investigation of stochastic effects this property will be very 
useful. 

The universal form of g su b 1 is drawn in Figure 3. 

We see that the maximum of subintegral expression is at p = 3 and 
at p = 1.4 = pi, the half of the maximal height is attained. This value 
can be regarded as some characteristic boundary. The other boundary is 
p = p c = 5.8. 

4 First iteration with shifted maximum 

The value c = 0.189 corresponding to the maximum at x = was established 
after explicit numerical solution of One can pose a question what value 
of c corresponds to the maximum of the first iteration at x = 0. Certainly, 
this procedure will lead to c = 1/6 and to a shifted first iteration 

ft sh ~ exp(x - exp(x)) 

This iteration coincides with precise solution even better than the previous 
first iteration. It can be seen from Figure 4. There are two curves: the curve 
with sharp slope corresponds to the shifted first iteration and the curve with 
smooth slope corresponds to the precise solution. 
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5 Monodisperse approximation 

Let us calculate expression for g according to (JTJ). We shall use a monodis- 
perse approximation. We get 

g = c f (z-xf exp(a; - g(x)) w cNz z max 

J —oo 

where z max is the coordinate of a monodisperse peak and iV is the number 
of droplets. For N one can propose 

N = exp(x — g)dx 

J —oo 

For z max since the maximum occurs at z = —3 one can get 

Zmax Z -\- 2) 

This forms the first variant of monodisperse approximation. It is not very 
precise because we have attributed the value z + 3 to small droplets which 
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Figure 4: The form of precise solution and the shifted first iteration 



are in a big quantity when we are near maximum of supersaturation. So, we 
have to correct our approximation. 

The region of integration to calculate iV can be cut off at pb- Then the 
renormalized number of droplets formed until maximum of supersaturation 
is 

N(z = 0) = / exp(x — g)dx 

J —oo 

For p > pb one can assume that g is negligible. This leads to the possibility 
to use the ideal supersaturation instead of the real one. Then one can come 
to 

roo 

N ml exp(—p)dp = exp(— p b ) = exp(— 1.4) 
Jpb 

So, at every moment the number of droplets is exp(— 1.4) exp(z) and the 
coordinate of monodisperse spectrum is x = 3. Then for / ~ exp(x — g) one 
can get 

/ ~ fmono i = exp(x - 3 3 cexp(-1.4) exp(x)) 

Figure 5 illustrates this approximation, here the upper curve is the precise 
solution and the lower curve is the current approximation. 
The precision of approximation is satisfactory. 



Figure 5: The form of monodisperse approximation /, 



This was only one variant of monodisperse approximation. Now we shall 
present another variant. Note that at every moment one can say that the 
droplet corresponding to the maximum of subintegral equation (in the first 
iteration) has a size p three units greater than a current size. This leads to 

/ ~ fmono 2 = exp(x - (x + 3) 3 CXp(-1.4)c) 

The form of spectrum is drawn in Figure 6. The smooth hill corresponds 
to the approximation, the sharp hill corresponds to the precise numerical 
solution. 

The precision of the last approximation is satisfactory. 

6 Monodisperse approximation at the maxi- 
mum of supersaturation 

Earlier we constructed monodisperse approximations for every moment of 
time. Now we shall take into account that the main quantity of droplets will 
be formed near the maximum of supersaturation, i.e. near z — 0. So, here 



Figure 6: The form of monodisperse approximation f mono 2 

we can put the position of monodisperse spectrum as z = —3. Then one can 
state that the maximum of the model spectrum 

/ ~ exp(x - A(x + 3) 3 ) 

with parameter A has maximum at z — 0. This gives A = 1/27. Moreover 
one can move back from the shifted spectrum to the unshifted one by the 
substitution of 6c/27 instead of 1/27. Then one can come to 

/ ~ fmono 3 = exp(x - — (x + 3) 3 ) 

This spectrum together with the precise solution is drawn in Figure 7. Here 
the sharp curve corresponds to the precise solution and the smooth curve 
corresponds to the approximation. It is seen that they are practically similar. 

When A = 1/27 one can come to 

/ ~ fmono 4 exp(x - ^(x + 3) 3 ) 

This spectrum together with the precise solution is drawn in Figure 8. Here 
the lower curve corresponds to precise solution and the upper curve - to 
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Figure 7: The form of monodisperse approximation /, 



approximation. It is seen that they are also practically similar, especially 
their forms. 

Instead of N ~ exp(— 1.4) exp(x) we put N ~ 6c/27 or N ~ 1/27 and 
conserve the fixed coordinate x — 3. 

The same can be done with approximations of the type f mon o i and f mon o 2 
Then we come to 

6c 

fmono 5 = exp(x - — exp(x)3 3 ) 

and 1 

fmono 6 = exp(x - — exp(x)3 3 ) 

We see that they coincide with f mono 3 and f m ono a- They are drawn in 
Figure 7 and Figure 8. So, we see that we can do the mentioned substitution 
at different steps of derivation. 



7 Two stage scheme 

The next approach to construct some approximations will be based on the 
property which results in applicability of the first iteration. What is the 
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Figure 8: The form of monodisperse approximation f mono 4 

physical reason of applicability of the first iteration instead of the precise 
solution? The physical reason is the following: the embryos which are the 
main consumers of vapor were formed at initial moments of nucleation period 
when the supersaturation was near the ideal value. This allows us to extract 
the initial stage of nucleation period z < — 1 and the second stage which is 
the whole nucleation period except the initial stage. The second stage is the 
stage when the main quantity of embryos was formed. 

Really g{— 1) < cexp(— 1) is small in comparison with characteristic value 
of g during the second stage g ~ c. So at z < — 1 the supersaturation ( is 
close to $. 

The influence of embryos appeared at initial stage on the evolution at the 
second stage is given through 



-1 



x 1 exp(x — cg(x))dx i — 0,1,2,3 



—00 



Direct calculation on a base of the first iteration leads to 



a = exp(-l) 



ai = —2 exp(— 1) 
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a 2 = 5 exp(— 1) 
a 3 = — 16exp( — 1) 

Then 

g ps Pin = c(x 3 a — 3x 2 ai + 3xa 2 — a 3 ) 

and the spectrum / ~ exp(x — g) — exp(a; — c(x 3 a — 3x 2 ai + 3xa 2 — a 3 )) is 
drawn in Figure 9 together with the precise solution which has more sharp 
form. 
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Figure 9: The form of approximation of two stages together with precise 
solution. 

Here we forget that this expression takes place only for z > —1. We 
spread it for z < — 1 because here the spectrum isn't too high. 

However, the asymptotic behaviour at z — > — oo will be qualitatively 
wrong. But as for the form of spectrum in a really significant region we see 
that the coincidence of two curves is practically ideal. 

8 Balance property 

The model forms of the size spectrum obtained above have to be self con- 
sistent. The self consistency can be treated at least in two senses. The 
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first is the weak dependence of the total number of droplets on the number 
of droplets in the monodisperse peak. This question is separately analyzed 
directly in investigation of stochastic effects. Here we shall see the weak de- 
pendence of the total number of droplets on the position of the monodisperse 
peak or on the position of the boundary between two stages of the nucleation 
period. 

We shall analyze two models - the model of the monodisperse peak and 
the model appeared in the two stage scheme. These models are the most 
precise ones. 

We shall start with the model of the monodisperse peak. At first we shall 
generalize it. Let the monodisperse peak be formed at z — —p. Then instead 
of A = 1/27 we have A = l/(3p 2 ). For the size spectrum we get 



The total number of droplets we get the following expression 

N t ot = N in + exp(-p) 

where N in is the number of droplets appeared after z = —p. The last value 
can be calculated as 



In the renormalization adopted here the total number calculated on the 
base of precise solution is precisely equal to 1.00. The numbers N in and N tot 
are drawn in Figure 10 as functions of p. It is clear that the dependence of 
Ntot on p is very weak one. 

The same we shall do for the model of the two stage approach. We shall 
state that before the boundary z = —p the spectrum is the ideal one 



Fot z greater than — p the evolution is governed only by the substance accu- 
mulated in the droplets appeared before z = —p. This gives 



/ ~ exp(x 



1 



3p 




f ~ exp(x) 



/ ~ exp(x — c(x 3 a — 3x 2 ai + 3xa 2 — 03)) 



where universal momentums are given by 



a = exp(-p) 
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at = —pexp(-p) 
a 2 = p 2 exp(— p) + 2pexp(— p) + 2exp(— p) 
03 = — p 3 exp(— p) — 3p 2 exp( — p) — 6pexp(— p) — 6exp(— p) 
The total numbers of droplets N to t is equal to 

N tot = exp(-p) + N in 

where iVj„ is the number of droplets appeared after z = —p. The last value 
is given by 

/°° 
exp(x — c(x 3 ao — 3x 2 ai + 3xa 2 — a^))dx 
-p 

The numbers N in and N tot are drawn in Figure 11 as functions of param- 
eter p of the boundary z = —p. One can see that the dependence of N tot on p 
is very weak one. Moreover the number of droplets (which is always greater 
than the precise value 1.00) has minimum at p ~ 1. Namely this value is the 
optimal one. We see that the value p = used in [3] is very far from the 
optimal one. 
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One can also note that this property of the weak dependence in deeply 
associated with constructions in the base of the modified gaussian method 
considered in 0. 
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Figure 11: The values N tot and N in as functions of the position of the 
boundary p. 
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